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The HII regions around quasars and galaxies at redshifts beyond the epoch of 
reionisation will provide prime targets for upcoming 21cm campaigns using a new 
generation of low-frequency radio observatories. Here we show that the boundaries of 
these HII regions will not be sharp. Rather, the clustering of sources near massive 
galaxies results in a neutral fraction that rises gradually towards large radii from an 
interior value near zero. A neutral fraction corresponding to the global background 
value is typically reached at a distance of 2-5 times the radius of the HII region around 
the central massive galaxy. 
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1 INTRODUCTION 

The identification of a Gunn-Peterson trough and ionised 
bubbles in the spectra of quasars at redshifts of z ~ 6.3-6.4 
(Fan et al. 2004) hints that the reionisation of cosmic hy- 
drogen was completed only around that time, about a bil- 
lion years after the big bang (White et al. 2003; Wyithe & 
Loeb 2004). In the near future, the best probe of the epoch 
beyond reionisation will be provided by a new generation 
of low-frequency radio arrays that will image the redshifted 
21cm emission from the neutral hydrogen in the intergalactic 
medium (IGM) 1 . The first target for these observations will 
be provided by the cosmological HII regions generated by 
massive high redshift galaxies and quasars. These cosmolog- 
ical HII regions are generally modeled with sharp boundaries 
within an IGM having a global neutral hydrogen fraction. 
However clustering of small galaxies around massive systems 
could lead to enhanced ionisation in the IGM immediately 
surrounding HII regions. In this Letter we estimate this ef- 
fect. Our approach is to use our earlier model for HII regions 
around clustered sources (Wyithe & Loeb 2005; hereafter 
Paper I). Throughout the paper we adopt the set of cos- 
mological parameters determined by WMAP (Spergel et al. 
2006) for a flat ACDM universe. 



2 THE SIZE OF IONISED REGIONS AROUND 
ISOLATED GALAXIES 

We begin by summarising the features of our model for HII 
regions around clustered galaxies. Only a brief description 
is given here, and we direct the reader to Paper I for further 
details. There are two contributions to the ionising luminos- 
ity of a galaxy, which we describe in turn: massive stars and 
a central supermassive black-hole. 

In a neutral IGM the characteristic physical size of a 
spherical HII region generated by stars in a galaxy with a 
dark matter halo of mass Mhaio is given by (Loeb et al. 2005) 
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where /* is the fraction of baryons in the halo that are turned 
into stars, and / csc is the fraction of ionising photons that 
escape from the galaxy into the surrounding IGM. The num- 
ber of ionising photons per baryon incorporated into Pop-II 
stars is N* ~ 4000 (Bromm, Kudritzki & Loeb 2001), and 
iVreion denotes the number of photons required per baryon 
for the reionisation of the HII region. The value of N le i on 
depends on the number of recombinations and hence on the 
dumpiness of the IGM. The halo mass Mhaio is a function of 
its virial velocity dispersion a, which we set equal to v c /^/2, 
where v c is the virial circular velocity [see Eqs. (22)-(25) in 
Barkana & Loeb (2001)]. 

Similarly, the physical size of the HII region generated 
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by an accreting supermassive black-hole of mass Mbh is 
(White et aL, 2003) 
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where iV 7i 5 is the production rate of ionising photons dur- 
ing the luminous quasar phase for a black hole of mass 
Mbh = 1O 5 M0, and /d yn = 0.035 is the fraction of the 
halo dynamical time during which the quasar shines at its 
peak luminosity. The black hole mass is related to the halo 
velocity dispersion through Mbh = 10 5 (er/54 km s -1 ) 5 Mq 
(Wyithe & Loeb 2003). 



3 THE STRUCTURE OF HII REGIONS 
AROUND CLUSTERED GALAXIES 

Prior to reionisation, isolated galaxies generate their own 
Stromgren (HII) spheres with characteristic sizes given by 
equations Q and @. However, massive galaxies tend to 
form in high density environments within which many 
smaller galaxies are clustered and contribute to the same 
HII region. 

Fixing the center of coordinates to be at the center of 
largest local galaxy, we may use equations Q and (2) to find 
the physical size i? s of the HII region surrounding a galaxy 
of mass Mhaio- We obtain 
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where the radii of the HII regions around the central halo 
of mass Mhaio and a galaxy of mass M, respectively, are 
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tion (pj determines R s self-consistently with the mean neu- 
tral fraction into which the HII region expands, and accounts 
for the higher IGM density due to infall around a massive 
halo through the nonlinear density contrast Ar(J?, Mhaio). 
The number of galaxies surrounding the central halo is 
found from the Press-Schechter (1974) mass function of ha- 
los dn(R, Mhaio) / dM with the modification of Sheth & Tor- 
man (2002), expressed in units of (physical Mpc) _3 Mg . 

In the vicinity of a massive halo dn/dM must be mod- 
ified relative to the background universe (see Paper I). The 
critical value of linear overdensity at which the actual col- 
lapse of a spherical shell occurs (<5 C rit = 1-69 for fi m = 1) 
depends weakly on the background matter density (Bryan & 
Norman 1998). The dominant contribution to the modifica- 
tion of the halo mass function originates from the accelerated 
formation of halos in overdense regions. 

In the background universe a spherical top-hat collapses 
at redshift z if its overdensity extrapolated (linearly) to the 
present day reaches 8 C iit(z) = 1.69[D(z)/D(0)] -1 , where the 
function D(z) is the growth factor. In the overdense region 
surrounding a massive galaxy the growth factor is modified 
relative to the background universe because of the enhanced 
matter density. Each spherical shell surrounding the central 
halo behaves as if it belongs to a universe with a modified 
value of the density parameter Q m given by its radial ex- 
pansion rate and the average matter density at the radius 



of the shell. We use the corresponding value of the growth 
factor for the value of Q m associated with each shell. 

Within the HII region the neutral fraction is close to 
zero (zhi ~ 0). Beyond the HII region, the clustering of 
galaxies may still yield ionisation at a level that is larger 
than the global average. The structure of ionisation in the 
IGM prior to reionisation is thought to be like that of a 
sponge, the neutral IGM being punctuated with HII re- 
gions (each with a sharp boundary) of different sizes. An 
observation at infinite resolution would be able to resolve 
all these HII regions. However a real instrument will ob- 
serve this sponge-like structure at finite resolution. Beyond 
the radius of the HII region we therefore define the local 
neutral fraction to be the mass fraction of neutral hydrogen 
averaged within a region that is smaller than the size of the 
main HII region, but large enough to contain many HII re- 
gions around small galaxies. This is the appropriate measure 
to consider for observations at finite spatial resolution (e.g. 
redshifted 21cm observations). 

At a radius R from the central massive galaxy, the neu- 
tral hydrogen fraction is 
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At large R, dn (R, Mhaio) / dM approaches the average mass- 
function in the Universe, Ar approaches unity, and equa- 
tion Q reduces to the average neutral fraction in the IGM 
(see equation^. Thus, away from the central HII region the 
above definition of neutral fraction is equivalent to the mean 
neutral fraction of the IGM. 



4 MODEL PARAMETERS AND THE 
REIONISATION HISTORY 

In order to compute the size and profile of an HII region, 
we need to specify the star formation efficiency and its de- 
pendence on halo mass, as well as the escape fraction of 
ionising photons from galaxies. As in Paper I, we define four 
models which bracket the expected range of star formation 
efficiencies. We then fix the free parameters in these models 
by requiring that the resulting reionisation history match 
observational constraints. 



4.1 Models for the star formation efficiency 

First we describe the four different models for the star 
formation efficiency that are examined in this paper. We 
adopt feedback-regulated prescriptions (Dekel & Woo 2002; 
Wyithe & Loeb 2003) in which the star formation efficiency 
in a galaxy depends on its velocity dispersion a, with f* = 
for a < <7 m m, /* = /*, cr it (cr/(j cr it) 2 for a min < a < a CX it, and 
/* = /*, C rit for a > (J cl it. We consider four cases in this pa- 
per: (A) a min = <7 cr it = 10 km s -1 ; (B) cr min = 10 km s _1 , 
Ocrit = 120 km s _1 ; (C) cr m m = <7 C rit = 50 km s~ x ; and 
(D) <T m i n = 50 km s _1 , CT cr it = 120 km s~ x . The values of 
o"min = 50 km s _1 , and <r m i n = 10 km s _1 , correspond to 
the threshold for infall of gas from a photo-ionised IGM and 
the threshold for atomic hydrogen cooling, respectively. The 
value of (Tcrit = 120 km s _1 is calibrated empirically based 
on the local sample of galaxies (Kauffmann et al. 2003; Dekel 
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Figure 1. The mean neutral fraction of the IGM, xhi, as a function of redshift z for the four models of star formation A-D (see text 
for details). 



& Woo 2002). Case A is a fiducial model where all galax- 
ies have the same star formation efficiency, and form down 
to masses corresponding to the hydrogen cooling threshold. 
This model describes the situation at high redshift where 
the IGM was not pre-heated, and where feedback does not 
operate in small galaxies as at low redshift because the dy- 
namical time within galaxies is comparable to the lifetime 
of massive stars. Cases B-D introduce modifications due to 
galaxy formation in a heated IGM and to stellar feedback. 

4.2 The reionisation history 

As described in Paper I, we may constrain the free param- 
eters in each of the four models defined above by requiring 
that hydrogen be reionised by z = 6 (Fan et al. 2004), and 
helium be doubly reionised by z = 3.5 (Theuns et al. 2002). 
The global ionisation fraction of intergalactic hydrogen is 

,-, f°° r dn 4-7T 

,mi = J Mmin 1m Y (5) 

Reionisation histories based on equation @ for the different 
star formation prescriptions A-D which were described in 
the previous section are plotted in Figure^ Note that these 
histories (and subsequent results) are independent of N Te i on 
which is degenerate with other parameters (e.g. /star)- 

The total optical depth for electron scattering is r os = 
0.06 in case A and ~ 0.05 in cases B-D. All four models are 
consistent with the latest measurement of r cs = 0.09 ± 0.03 
from WMAP (Spergel et al. 2006). The error bars on the 
WMAP data do not allow us to distinguish between the 
models. 



5 RESULTS 

We now examine our results for the effect of clustered 
sources on the neutral fraction outside the HII regions sur- 
rounding massive galaxies (equation |1J. 

5.1 The neutral fraction outside quasar HII 
regions 

In Figure [5] we plot the neutral fraction xm as a func- 
tion of radius R around a central quasar. The four panels 
in Figure [5] refer to cases A-D, and the solid and dashed 
lines correspond to velocity dispersions of 350 km s" 1 and 



200 km s~ for the central galaxy, respectively. In each case 
the four sets of curves correspond to redshifts of z = 7.5, 
6.9, 6.3 and 6.1 (top to bottom). These calculations of 
xm(R) are self-consistent with the average neutral fraction 
of the IGM, which is obtained from the reionisation histories 
(equation |5J|. The figure shows that the edges of HII regions 
are not expected to be sharp, particularly when the neutral 
fraction is of order unity. Indeed, in such cases the neutral 
fraction rises gradually beyond Rs, and reaches the global 
value only at radii that are 2-5 times R s . 

Figure [5] shows that the detailed profile of the neu- 
tral fraction with radius depends on the value of the global 
neutral fraction. In all cases the global neutral fraction is 
reached at a radius ~ 8Mpc. However galaxies produce 
larger HII regions within a more ionised IGM. As a result, 
the boundary of the HII region is spread over a larger range 
of radii (relative to R s ) in cases where the global neutral frac- 
tion is higher. The profile is also sensitive to the assumed 
case of star formation. In particular, consider star-formation 
scenarios where there is only a small contribution to the 
ionising photon production from galaxies with a < 50km/s. 
The boundary is less extended in such cases, particularly 
where the HII region forms within an IGM that is already 
significantly reionised. On the other hand, for large global 
neutral fractions the curves corresponding to the four mod- 
els are similar. This is simply because the parameters in 
these models have been tuned to provide hydrogen reionisa- 
tion at z ~ 6. By pinning the models to the same redshift of 
reionisation, we draw the robust conclusion that HII regions 
around high redshift quasars within a substantially neutral 
IGM should not have sharp edges, but rather the neutral 
fraction in the IGM will gradually reach the global value. 

We note that the effect on the ionisation fraction is 
significant relative to the background level of ionisation on 
scales as large as ~30 comoving Mpc. Although the over- 
density on this scale is mild, ~ 10%, it nevertheless biases 
the abundances of galaxies, particularly those on the expo- 
nential tail of the mass function, to form earlier (Barkana 
& Loeb 2004). This early formation produces an ionisa- 
tion level that is significantly larger than that of the back- 
ground. The magnitude of the effect we calculate agrees with 
expectations for the clustering of galaxies (Scannapieco & 
Barkana 2002). For example, given a halo mass of ~ 1O 9 M 
(that is a factor ~ 10 3 smaller than a central galaxy with 
a — 350km/s at z = 7.5), we find the cross-correlation 
length with the central halo to be ~ 1.3 physical Mpc. This 
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Figure 2. Neutral fraction as a function of radius around a quasar at z > 6. The four panels show results for cases A D. In each case 
the solid and dashed curves show central velocity dispersions of a = 350 km s — 1 and a = 200 km s _1 respectively. The four sets of 
curves from top to bottom correspond to redshifts of z = 7.5, 6.9, 6.3 and 6.1. 



length represents density enhancements of order unity and is 
naturally smaller than the ~ 3-5Mpc where enhanced ioni- 
sation is seen outside the HII region. However we find a 50% 
enhancement in the density of halos at ~ 2Mpc (physical 
scale) relative to a random distribution, and a 10% enhance- 
ment at ~ 4Mpc. For a central galaxy with a — 200km/s, we 
find a correlation length of ~ 0.8Mpc, with an enhancement 
of 50% at a distance of ~ 1.3Mpc and of 10% at ~ 3Mpc. 
These enhancements can naturally account for the 10-30% 
increase in the ionisation level at these radii. 



5.2 The neutral fraction outside galaxy HII 
regions 

In Paper I we demonstrated that quasars do not generally 
dominate the contribution to ionisation within an HII region 
around clustered galaxies. Therefore even in the absence of 
quasar activity we would expect a large HII region surround- 
ing a massive galaxy. Clustering of galaxies around a central 
massive galaxy should also lead to a gradual boundary of the 
HII region as is the case for a central quasar. 

We therefore repeat the calculations of the previous sec- 
tion in the case where quasar activity is ignored (both in 
the central galaxy, and surrounding galaxies), by replacing 
£>Lio = ^haio.q + £>Lio,* and D M = D M , q + D M>ir , with 
Aiaio = Aiaio,* an d Dm = Dm,* in equations The 
results are shown in Figure^] We find profiles that are quali- 
tatively consistent with those shown in Figure[5] The quanti- 
tative difference is that HII regions around massive galaxies 
are smaller in the absence of quasars. As before, the neutral 
fraction rises gradually beyond _R S , and reaches the global 
value only at radii that are 2-5 times i? s . 



6 SUMMARY 

Our calculations suggest that observations of cosmological 
HII regions will not see sharp boundaries, but rather a grad- 
ual rise in the average local neutral fraction from the radius 
of the HII region out to a distance 2-5 times greater. This 
gradual rise in neutral fraction is produced by the clustering 
of sources around the massive systems at the centers of cos- 
mological HII regions. The gradual rise is seen in our models 
for HII regions surrounding high redshift quasars, as well as 
massive galaxies in the absence of quasar activity. The grad- 
ual nature of the HII region boundary is most prominent in 
star formation scenarios that do not include strong feedback 
in low mass galaxies. 

Our Letter only considered the softening of the HII re- 
gion boundary by galaxy clustering. Additional softening of 
the boundary could be produced near quasars that emit sig- 
nificant energy in X-ray photons, whose mean-free-path is 
much longer than UV photons. For an unobscured quasar 
the spectrum averaged mean-free-path is much smaller than 
the size of the HII region. As a result, the high energy pho- 
tons smooth the edge of the HII region over a length much 
smaller than its radius. Along directions where the UV emis- 
sion is obscured but the X-rays are not, the shape of the HII 
boundary might be dominated by the quasar spectrum and 
not the clustering of galaxies as discussed here. 

We have assumed spherical symmetry in our model. 
However, in reality massive galaxies and quasars will lie at 
the intersections of filaments along which mass is tunneled 
into them. The shape of the HII region and the level to 
which the boundary is gradual, will vary along different di- 
rections. Numerical simulations are needed to explore this 
anisotropy. Nevertheless, the first observations of HII regions 
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Figure 3. Neutral fraction as a function of radius around a massive galaxy at z > 6 in the absence of quasar activity. The four panels show 
results for cases A D. In each case the solid and dashed curves show central velocity dispersions of a = 350 km s" 1 and a = 200 km s — 1 
respectively. The four sets of curves from top to bottom correspond to redshifts of z = 7.5, 6.9, 6.3 and 6.1. 



at redshifted 21cm wavelengths will be made at a resolution 
comparable to the size of the HII region in both transverse 
and line-of-sight directions (Wyithe, Barnes & Loeb 2005). 
This coarse resolution will tend to smooth the neutral dis- 
tribution in spherical shells as assumed in our analysis. 

The first generation of low-frequency arrays will likely 
detect HII regions via matched filter techniques. In the case 
of HII regions also observed through Lya absorption, the 
clustering of sources results in a "characteristic size" of an 
HII region (as measured assuming a sharp boundary) that 
is 1.5 to 2 times greater than the value measured via Lya 
absorption. Searches for HII regions using redshifted 21cm 
observations around known quasars should therefore use a 
filter that accounts for the gradual profile of the HII region 
boundary. 

The detection of the distribution of HII region sizes, to- 
gether with the power-spectrum of 21cm emission will pro- 
vide a detailed picture of the topology of the reionisation 
epoch (Furlanetto et al. 2004). We have demonstrated that 
there is a need to define the size of an HII region in a reso- 
lution dependent way. In particular, finite resolution obser- 
vations (for example using low-frequency radio telescopes 
to detect redshifted 21cm radiation) will not detect sharp 
boundaries around large HII regions. Rather, confusion of 
many small, clustered and unresolved HII regions, will blur 
the observed boundary of the HII regions around massive 
galaxies and quasars, and potentially increase their inferred 
HII volumes by an order of magnitude. 
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